Full counting statistics of Luttinger liquid conductor 
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Non-equilibrium bosonization technique is used to study current fluctuations of interacting elec- 
trons in a single-channel quantum wire representing a Luttinger liquid (LL) conductor. An exact 
expression for the full counting statistics of the transmitted charge is derived. It is given by Fred- 
holm determinant of the counting operator with a time dependent scattering phase. The result has a 
form of counting statistics of non-interacting particles with fractional charges, induced by scattering 
off the boundaries between the LL wire and the non-interacting leads. 

PACS numbers: 73.23.-b, 73.40.Gk, 73.50.Td 
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Fluctuations are among the most fundamental con- 
cepts arising in statistical physics. In recent years, non- 
equilibrium noise was measured in a variety of electronic 
systems, such as quantum point contacts diffusive 
mesoscopic conductors 01, and fractional quantum Hall 
edges [al; see Ref. i for review. With the decrease of a 
sample size, characterization of current fluctuations by 
the second moment only becomes insufficient. This has 
triggered recent interest in higher-order correlation func- 
tions of current statistics. The third cumulant of the 
noise was measured in recent experiments A more 
complete characterization of current fluctuations is the 
full counting statistics (FCS), introduced by Levitov et 
al. 0. In its full extent this fascinating theoretical ap- 
proach yields information about all moments of the num- 
ber of electrons transferred (over a given time interval) 
through a given terminal in a multi-terminal system, cur- 
rent cross-correlations and entanglement, and large cur- 
rent fluctuations. 

For non-interacting systems, the FCS of non- 
equilibrium fluctuations is well understood by now within 
several complementary approaches, including the Fred- 
holm determinant formalism 01 , the cr-model field- 
theoretical description as well as the kinetic theory 
of fluctuations and related methods [9|. 

Much less remains known concerning fluctuations in 
interacting systems. This problem is of particular in- 
terest at low dimensions, where interaction affects the 
nature of the system at hand in an essential way. This 
is, in particular, the case for FCS of current through a 
quantum dot in the Coulomb blockade regime The 
interaction also affects dramatically the physics of one- 
dimcnsional (ID) fermionic systems where a strongly cor- 
related state — Luttinger liquid (LL) — is formed[ll|. Ex- 
perimental realizations of LL include carbon nanotubes, 
semiconductor, metallic and polymer nanowires, as well 
as quantum Hall edges. Recent experiments studied the 
shot noise in carbon nanotubes (l2l |. Non-equilibrium 



physics of carbon nanotubes and quantum Hall edges 
has been explored through tunneling spectroscopy [13| 
and Mach-Zehnder-iterferometry |14] respectively. The 
LL description is also relevant to interacting ID bosonic 
systems, cf. e.g. recent experiments on ultracold atomic 
gases that probe statistical properties of inteference con- 
trast, thus obtaining information on the full distribution 
of quantum noise [l5|. 

Previous theoretical work on current fluctuations in LL 
mainly focussed on the second moment (shot noise) [l6| . 
The most intriguing observation was related to manifes- 
tations of fractional charges in shot noise. However, full 
understanding of the nature of charge transfer processes 
(and, in particular, of charge fractionalization) requires 
the analysis of the FCS. The latter has been studied in a 
biased LL with an impurity [l3l ■ While an analytical so- 
lution via the thermodynamic Bethe ansatz can be found 
17| . it is in general implicit and very cumbersome. 

In the present work we study the statistical proper- 
ties of current fluctuations in non-equilibrium LL con- 
ductor with arbitrary energy distribution of right- and 
left-movers. Employing the recently developed non- 
equilibrium bosonization technique [18|, we find an ex- 
act solution of this problem. We show that the full dis- 
tribution of the current noise is closely related to the 
phenomenon of charge fractionalization. Specifically, we 
demonstrate that the FCS of LL reduces to the one of a 
non-interacting model with fractional charges calculated 
below. 

We consider a LL conductor geometry schematically 
shown in Fig.[T] Our goal is to solve FCS problem in the 
presence of electron interaction. We assume that the spa- 
tially dependent interaction strength g{x) vanishes out- 
side the interval [— L/2, L/2]; t his way of modeling leads 
was introduced in Refs. Il9l42ll to study the conductance 
of LL wire. To describe the FCS of charge transfer, we 
define a probability distribution function pr (n) for n elec- 
trons to pass through the cross-section during the time 
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FIG. 1: Schematic view of noise detection in LL conductor. 
The positions of the charge detector inside (A) and outside 
(B) of interacting region are indicated. The a;-dependent LL 
parameter K{x) for smooth (sohd hne) and sharp (dashed 
hne) boundaries is also shown. 



interval t. It is convenient to introduce a generating fuc- 

tionA«(A) = Er=-ooPrHe"'^ 

Let us first discuss the non-interacting case. In that 
case k(A) has been calculated in Ref. by means of Lan- 
dauer scattering-state approach. For an ideal quantum 
wire (with no scattering inside the wire) the generating 
function of FCS is given by 

k{X) ^ AR[SR{t)]AL[SL{t)] . (1) 

Here t] — R,L labels right-and left-movers and 

A^[<5,(t)] =det[l + (e-*'" - IK] (2) 

is a Fredholm determinant (of Toeplitz type) of the 
counting operator built of fermionic distribution function 
n^(e) and a time-dependent scattering phase S,j{t), with 
e and t understood as canonically conjugate variables. In 
Eq. ([T]) and below we assume normalization of the deter- 
minant to its value for equilibrium, zero-temperature 
distribution. The phase Srj(t) is given by 



Sjjit) = Xr]Wr{t, 0) , 



(3) 



where we have defined a window function Wr{t,t) — 
9{t — t) — 9{t — t — t). We use the convention that in for- 
mulas 7] is understood as = ±1 for right-/left-movers. 
In the long time limit the Fredholm determinant © can 
be easily evaluated analytically; in a more general situ- 
ation it can be efficiently studied numerically Be- 
low we show that for the interacting case, the generating 
function of the FCS obeys the form of Eq. ([T]). All elec- 
tron interaction effects are encoded in the time dependent 
scattering phases (5^(i). 

On the microscopic level the problem is described by 
the Keldysh action S = Sq [t/j] + Sac [V'] j written in terms 
of fermionic fields 2J] , 



^ Jc J 



Here = V'JV'r; are density fields, dj^ = dt + rjvdx, 
V is the Fermi velocity. The non-equilibrium bosoniza- 
tion approach allows us to reformulate this theory in 
terms of bosonic (density) fields. The interacting part 
of the action. See , is already expressed in terms of den- 
sity modes p^. Following rotation in Keldysh space 
See — J2riri' I dtdxg{x) Pripri' , whilc the free part of the 
action, where information concerning the state of the 
non-interacting fermionic system is encoded, reads 



(4) 



Here we have decomposed the bosonic variables into clas- 
sical and quantum components, p,p — {p+ ± p_)/a/2, 
where the indices -|- and — refer to the two branches of 
the Keldysh contour; HJJ is the advanced component of 
polarization operator, and Z^lx-q] is a partition function 
of free fermions moving in the field 

Xv = 'Pr, , n^(w, q) = T]q/2TT{r]vq - uj - iO) . (5) 

Expansion of Z^\xr\ in x^(i,a;) generates an infinite 
series, ilnZ^[xr,] = E„(-1)""^^XJJ'5^"/"; governed 
by irreducible fermionic density correlation functions, 
5,,„(ti,a;i; . . . ;t„,a;„) = ((p,,ip2 • ■ • Pr,n)), representing 
cumulants of quantum noise 6]. 

To find the generating function for the FCS of charge 
transferred through a cross-section x = X during a time 
interval r, one needs to calculate the correlation function 



«:(A,X) = (e'^Q(^'^)e-^^Q(''^^)) 



(6) 



where Q(t,X) — J^^ipnix) + pL{x))dx is the electric 
charge situated on the left side of the point X at the time 
t. We note a mathematical similarity between the prob- 
lem of FCS and that of tunneling spectroscopy The 
latter amounts to evaluation of a single-particle Green 
function. In both cases one needs to find a correlation 
function of exponentials of bosonic operators that are lin- 
ear combinations of the right and left density fields. A 
major difference stems from the fact that the operators 
gizAQ jj-^ ijgji fjQj-^^ain a sum of left and right densities 
in the exponential, whereas the fermionic operators tpri, 
-0^ contain only p,,. 

We now proceed with calculation of the generating 
function ([5|); technical details are outlined in the Sup- 
plementary Material [l^. We find that the FCS gener- 
ating function has the form of Eq. (HI), as in the case 
of free fermions. The time dependent scattering phases 
6,-j (t) are expressed through plasmon refiection and trans- 
mission coefficients r,j, at the left (I/H; rj = L) and 
right (II/III; rj — R) boundaries; r'^ + = 1. In gen- 
eral, multiple scattering off these boundaries will give 
rise to an infinite sequence of scattering phase pulses 
[cf. Eq. If the boundary is smooth (on the scale 

of the plasmon wave length), there is no plasmon reflec- 
tion, ~ 0. In the opposite limit of a sharp boundary 
we have = (1 - K)/{\ + K). 
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FIG. 2: Phases 5^ entering Eq. for the measuring device 
located at X = 0, for sharp (a,b) and adiabatic (c,d) bound- 
aries {K = 1/3). 



In contrast with the free fermion case, the statistics 
of current fluctuations in the LL conductor depend on 
the position of the measuring device (see FigH]). Let us 
first consider the case where the current fluctuations are 
measured in the middle of an interacting region {X = 0, 
position A in FiglT]). We thus find [23] that time depen- 
dent scattering phases consist of a sequence of pulses 



(7) 



with partial phase shifts 

<5r,,2n = r]Xt-,^VKrl^r'\ = ?/Ae* 2„ , (8) 

'5„,2n+i = ?7At-^Vifr;+V!!,, = vK,2n+i ■ (9) 

The beginning of the n-th pulse occurs at time t„ = 
{n + 1/2 — l/2K)L/u. The sequence of pulses 5r^{t) is 
shown in Fig. [2] for the limiting cases of sharp and smooth 
boundaries between the interacting and non-interacting 
regions. 

In the limit of very low frequencies (r ^ L/u) all 
pulses (O overlap and the scattering phase can be ap- 
proximated by 5n{i) — Wr(t, 0) Yl^=o ^v,n- In this limit of 
long wavelength the interacting/non-interacting bound- 
aries appear to be sharp. This yields X]^o ^v,n — V-^- 
As expected l^- 21 1 the effects of interaction in this limit 
disappear and one is back to the FCS of free fermionsQ. 

In the opposite case (r ^ L/u ) the interference be- 
tween different plasmon pulses may be neglected and the 
result splits into a product 



-n^eln . (10) 



with each factor representing a contribution of a sin- 
gle phase pulse <5,,^„(t) = 5^„i(;T-(t,0). Remarkably, 
plasmon scattering gives rise to charge fractionalization, 
which splits the scattering phase into an infinite series 
of pulses. On a technical level, the Fredholm determi- 
nant of the counting operator is now replaced by an 
infinite product of determinants, each calculated with 
a corresponding scattering phase. The FCS of the LL 
then becomes a superposition of FCS of non-interacting 
electrons with fractional charges e* „ . For the case of 
smooth boundaries we get only one fractional charge. 



e* Q = \JK. In the opposite limit of sharp boundaries, 
we obtain the sequence of fractional charges of the form 
e;,„ = 2K(l-K)"/(l + K)"+i. 

Focussing on the regime where the counting interval 
r is long compared to the inverse energy scale of the 
distributions rirfit), we express A,,r(^) in the form 



In A^,(J) = r / — {ln[l (e"^^ - l)7i„(e)] + i8B{-t)}. 
J 27r 

(11) 

Substituting ([TT]) in Eqs. ([!]), pHI) and expanding in 5, 
one can find explicitly all moments of noise. 

Consider now the FCS outside the interacting region, 
e.g. at X < -L/2 (B in FigH]). The analysis follows 
closely the one outlined above, hence we present only the 
final result. The right and left scattering phases are 



5R{t) = A 



, X\ ( L + X 

Wr [t,— j - rLWr I <, 



.2^ n n (L + X 2(n+ 1)L 



6L{t) - -XtLtRY^^l^RWr U + L 



vu 



Current fluctuations measured in the non-interacting re- 
gion at finite frequency (r ^ L/u) differ from those at 
the interacting part of LL ([5]). Here too the FCS takes 
the form of superposition of non-interacting FCS's, but 
with a new set of fractional charges: 



= 1, 



tLtRvlrl 



'R,n+2 — l^L^L^R 



0,1,...). 



(12) 



ri=0 



In the universal ultra- low- frequency limit (r ^ L/u) the 
effects of interaction again disappear, and the result for 
the FCS of free fermions is recovered. 

To conclude we provide a brief overview of the physics 
discussed here. We have calculated the FCS in a LL 
conductor. The result is cast in terms of a Fredholm de- 
terminant of the counting operator. The boundaries be- 
tween the interacting and non-interacting regions of the 
wire give rise to (in general, multiple) plasmon scattering. 
This generates plasmon wavelets corresponding to frac- 
tionally charged "particles". As a result, for a counting 
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time interval r small compared with the plasmon time- 
of- flight through the interacting region, L/u, the FCS is 
a superposition of counting statistics of non-interacting 
particles with fractional charges e* Let us stress that 
the FCS contains more complete information about the 
system than the second moment (noise) . As an example, 
for a symmetric system (r^ = rn) the noise acquires the 
factor K due to interaction. By itself, this is not suffi- 
cient to make a conclusion about the value of fractional 
charges and the character of correlations between them. 

The fractionalization process that manifests itself 
in the FCS caculated above is solely due to the 
interacting/non-interacting boundary scattering [20| . 
Fractionalization due to tunneling into a LL [25j does 
not show up here: technically this is since the operator 
Q in Eq. ([6]) is determined by a sum of pL and pR. 

How is the notion of charge fractionalization compat- 
ible with charge quantization? The above analysis was 
based on the bosonization approach with the density 
slowly varying on the scale of the inverse Fermi wave 
vector kp^. This corresponds to the situation where the 
measurement procedure is smooth on the scale kp^. In 
this situation the above results are valid for any counting 
field A: our detector may count a fractional charge resid- 
ing in some volume of space. In the opposite limit, where 
the spatial resolution of our measurement is sharp on the 
scale kp^, the measured charge should be integer, i.e. the 
FCS should satisfy k{X) = k{X + 27r). In the bosoniza- 
tion formalism, charge quantization emerges by taking 
into account fast oscillatory contributions to the density 
Pti This will not affect the moments of FCS (which 
are obtained as derivatives of k(A) at A = 0) and, more 
generally, the form of k{X) in the range [— tt, tt]. Beyond 
this interval, the FCS will be continued periodically, in 
agreement with the charge quantization requirement. 
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US-Israel BSF, ISF, Minerva Foundation, DFG SPP 
1285, EU GEOMDISS, and Rosnauka 02.740.11.5072. 
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ONLINE SUPPLEMENTARY MATERIAL 

Here we present details of the calculation of the gen- 
erating function ([6]). Despite the fact that the bosonic 
action (|4|) contains terms of all orders, evaluation of the 
generating function can be performed exactly. The point 
is that the action ^ is linear with respect to the classi- 
cal component of the density field. Thus, integration 
with respect to can easily be performed. It gives rise 
to a 5-function condition for the quantum component 
of the density. 



ZTT ZTT 



dtPL - dx {v 



27r 



)PL 



9_ 

27r 



PR 



A 

^2^ 



Here jr{x,t) is a source term. 



jr{x,t)^S{x~X)[S{t-T)-6{ty 



J- (SI) 



(S2) 



It remains to solve Eqs. (|S1[) for p^, to calculate Xri ac- 
cording to Eq. ([S]) and then to substitute it into the sec- 
ond term of the bosonic action (|4]). This yields 



(S3) 



We thus find that the generating function of the FCS 
problem has the form of Eq. ([1]), as in the case of free 
fermions. The time dependent scattering phases Sr/it) 
are given by the integral of Xr; along a free-electron tra- 
jectory. 



/OO 
dixv{i-t,Vvi) 
-OO 



-27rA/2 rj lim 



dx pri{x,t) . (S4) 



The second line of Eq. ([S4|) expresses S,j{t) through the 
asymptotic behavior of pri{x,t) in the non-interacting 
parts of the wire (regions I and HI in FiglT]). 



To find the scattering phases, we now solve the set of 
linear differential equations (IS1|) . These equations can be 
conveniently combined into a single differential equation 
for the current operator J ~ v{pfj — p^) 



{lu'^ + dxU^{x)dx)J{Lu,x) =0, x^X, 



(S5) 



where u{x) — v/K{x) is the plasmon velocity and 
K{x) = [1 + g{x)/TTv]~^^^ is the LL parameter. The lat- 
ter is equal to 1 in the leads and takes a constant value 
K in the interacting region; the crossover between these 
values depends on the profile of the interaction. Inside 
the interacting region J{uj, x) is given by a superposition 
of the right and left moving plasmon waves: 



Jiuj, x) = w-i/2(ae*'"^ + 66-'"^) , 



uj/u. (S6) 



In the non-interacting regions, an analogous form holds, 
with the replacement u ^ v. The change of interaction 
strength constitutes a source of plasmon scattering and 
can be described by plasmons scattering matrices. 



tL 

rL 



-rL 

tL 



Se 



-rR tR 



(S7) 



The matrices Sl,r relate the amplitudes a, 6 at the 
left (I/II) and right (II/III) boundaries, respectively; 
''^ + = 1- general, multiple scattering off these 
boundaries will give rise to an infinite sequence of scatter- 
ing phase pulses (cf. Eq. [7]). If the boundary is smooth 
(on the scale of the plasmon wave length), there is no 
plasmon reflection, ~ 0. In the opposite limit of a 
sharp boundary we have = (1 — K) /(I + K). 



Using Eqs. ([85]), ([86]), and we obtain the 

phases Sri{t) as given in the main text of the paper. 



